Doping dependence of the vortex glass and sublimation transitions in the high-T c 
superconductor La 2 _ a .Sr x Cu04 as determined from macroscopic measurements 
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Magnetization and ac-susceptibility measurements are used to characterize the mixed phase 
of the high-temperature cuprate superconductor I^-^Sr^CuCU over a large range of doping 
(0.075 < x < 0.20). The first order vortex lattice phase transition line Hfot(T), the upper 
critical field H C 2(T) and the second peak H sp (T) have been investigated up to high magnetic fields 
(8 Tesla applied perpendicular to the CuC>2 planes). Our results reveal a strong doping dependence 
of the magnetic phase diagram, which can mainly be explained by the increasing anisotropy with 
underdoping. Within our interpretation, the first order vortex lattice phase transition is due to the 
sublimation (rather than melting) of the vortex lattice into a gas of pancake vortices, whereas the 
second peak is related to the transition to a more disordered vortex glass state. 

PACS numbers: 74.25.Ha, 74.25.Qt, 74.25.Op, 74.72.Dn 



I. INTRODUCTION 

Despite belonging to the family of the first high-T c 
superconductor (HTSC) to be discovered, the magnetic 
phase diagram of I^-^Sr^CuO^ (LSCO) has not been as 
intensively investigated as that of other cuprates such as 
YBaaCusCX, (YBCO) and Bi 2 Sr 2 CaCu 2 8+;E (BSCCO). 
The LSCO compound has a relative small value of T c 
(38.5 K at optimal doping), but is of high interest because 
it fills the gap between rather 3D systems such as YBCO 
and highly anisotropic 2D systems such as BSCCO. The 
anisotropy factor 7 2 can be defined as the ratio between 
the out-of-plane and the in-plane resistive components 
[.Pel Pah) measured in the normal stated. An additional 
advantage of LSCO is that 7 depends on the Sr content 
x and allows a study of the magnetic phase diagram over 
a wide range of anisotropy (200 < j 2 < 4000) which lies 
inbetween the values for YBCO (25 < j 2 < 100) and 
BSCCO (3000 < 7 2 < 30000). 

The magnetic phase diagram of HTSC cuprates is dom- 
inated by the mixed phase (the lower critical field 
H cl (0 K) is about 10~ 2 T whereas the upper critical 
field -ff c2 (0 K) is of the order of 10 2 T), where the mag- 
netic flux can penetrate into the sample in the form of 
quantized flux- lines (vortices). Due to the anisotropy 
and thermal fluctuations one observes a number of vortex 
phases, which have been the subject of extended experi- 
mental and theoretical research in the last two decades^. 
In LSCO one can distinguish between a first order transi- 
tion (FOT) line Hfot(T), which has been attributed to 
the meltingMai or sublimation^^ of the vortex lattice into 
a vortex fluid, and the irreversibility line Hi rr (T), where 
reversible magnetization and resistivity appeariSiS. An- 
other interesting feature is the so-called fishtail effect, 
that is an anomalous second peak in the magnetiza- 
tion loops. The origin of the second peak line H sp (T) 
is controversial, and has been attributed to mechanisms 



varying from dimensional crossover"!, collective pinning 8 , 
crossover between different pinning phases, crossover to 
a disordered vortex glass2*i2*ii, etc. 

Only recently the vortex lattice (VL) has been directly 
observed in overdoped LSCO by means of small an- 
gle neutron scattering (SANS), revealing a field- induced 
transition from hexagonal to square symmetrjiiSii^ and 
the vanishing of the VL signal at temperatures well below 
T c2 i^. In the underdoped regime of LSCO, on the other 
hand, a more disordered vortex glass has been observed 
by means of muon spin rotation (/iSR) experiments^. 
Interestingly, recent inelastic neutron scattering (INS) 
experiments indicate a possible interplay between the 
vortex and copper-spin degrees of freedom. In opti- 
mally doped LSCO, sub-gap spin excitations induced 
by a magnetic field of 7.5 Tesla have been observed at 
low-temperatures 16 . Moreover, the spin gap was found 
to close at the irreversibility temperature rather than 
at T C 2 16 ' 17 . In underdoped LSCO, field-induced static 
incommensurate magnetic peaks have been observed 18 , 
and it has been suggested that these field-induced mag- 
netic signals arise from antiferromagnetic order in the 
vortex cores and in the surrounding regionsi2i22i2i. En- 
hanced antiferromagnetic spin correlations in the vor- 
tex core region have been indeed observed in NMR 
experiments^ 2 *^. 

In order to understand these experiments performed in 
the presence of an external magnetic field, it is crucial to 
have a good knowledge of the rich and complicated mag- 
netic phase diagram of HTSC. We will present here a 
detailed study of the doping dependence of the magnetic 
phase diagram in LSCO single crystals from a macro- 
scopic point of view. 
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II. EXPERIMENTAL 

The magnetic phase diagram of LSCO has been in- 
vestigated by means of magnetization (M) and ac- 
susceptibility (x) measurements. We used a Quan- 
tum Design Physical Properties Measurements System 
(PPMS) up to fields of 8 T applied approximately per- 
pendicular to the CuOi planes. The angle between the 
field direction and the c-axis of the samples was always 
smaller than 10 degrees. This precision is good enough 
for the present study, since the critical lines (e.g. melt- 
ing line H m , upper critical field H C 2) are known to be 
only slightly affected by small angles (e.g. H m (Q) ~ 
H m (Q = 0)/cos(Q), H c2 (Q) ~ H c2 (B = 0)/cos(9))» 
Four high quality LSCO single crystals with different 
doping levels have been measured. Details of the sample 
growth can be found elsewhere2&. The samples are la- 
beled by the doping region (OD for overdoped and UD for 
underdoped) together with their T c , defined by x'(^c) = 
i%'(0 K). The width of the superconducting transition 
AT C has been determined by the 10%-90% criterion. OD- 
31K is a highly overdoped (x=0.20, T c =31.5 K, AT C =2.8 
K) 51 mg crystal. OD-36K is slightly overdoped (x=0.17, 
T c =36.2 K, AT C =1.5 K) and is a portion of the crystal 
used for our SANS and INS experimentsiSiii. While most 
of the measurements on the OD-36K sample have been 
performed on a 293 mg cylindrical crystal, for zero-field 
cooled magnetization measurements the crystal has been 
cut to a 84 mg plate-like shape with the c-axis parallel 
to the largest face, in order to reduce the diamagnetic 
signal. UD-29K is an underdoped (x=0.10, T c =29.2 K, 
AT C =1.3 K) 37 mg plate- like crystal with the c-axis par- 
allel to the largest face, which has been cut from a larger 
crystal used in /xSR experiments^. Finally, UD-19K is a 
highly underdoped (x=0.075, T c =19 K, AT C =3.8 K) 52 
mg plate-like crystal with the c-axis perpendicular to the 
largest face. 



III. RESULTS 

We start with the complex ac-susceptibility x—x'+^x" ■ 
The samples are placed in an external magnetic field 
Hext = Hdc + H ac ■ cos(uJact), with H ac =10 Oe and 
uj ac = 10 Hz, T< H dc <8 T. A set of field-cooled (FC) 
temperature scans x{T) for the four LSCO samples in 
different magnetic fields is shown in FigQ] with the real 
part x' an d imaginary part x" plotted separately. In all 
samples the peak in x" shifts toward low temperatures 
and sharpens with increasing magnetic field. However, 
the magnitude of the shift is strongly doping dependent: 
for UD-19K a magnetic field of 6 T is sufficient to shift 
the peak by 0.85-T c , whereas for OD-31K the shift caused 
by a field of 8 T is only 0.45-T c . The detailed field depen- 
dence will be discussed in Section llVl In Fig[2]a represen- 
tative curve x{T) measured at Hd c =3 T for UD-29K is 
plotted together with magnetization curves M(T). One 
can notice that there is no difference between the zero- 
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FIG. 1: Real and imaginary part of the ac-susceptibility x{T) 
for OD-31K, OD-36K, UD-29K and UD-19K measured at dif- 
ferent magnetic fields between T and 8 T. The peak in x' (T) 
rapidly shifts toward lower temperatures with increasing field. 



field cooled (ZFC) and the FC x{T) data, whereas FC 
and ZFC M(T) curves separate below the irreversibility 
temperature Tj rr . Slightly above Ti rr there is a jump in 
the magnetization, indicating the presence of a first order 
transition (FOT)i*2. Similar data have been obtained for 
the other samples and for other values of Hdc- The jump 
is more pronounced at high magnetic fields, and in UD- 
I9K only a broad anomaly could be observed (to note 
is that in this sample the loss peaks in x"(T) are very 
broad, as well). 

The experimental Ti rr is often obtained from the loss 
peak in x"(T), which is directly related to the maximum 
slope in x'(T)2£. However, in our case, Ti rr obtained 
by ac-susceptibility measurements is slightly higher than 
the "real" Ti rri and is concomitant to the jump in M(T) 
at Tfot- The irreversibility line and the FOT line are 
found to be close to each other in all the samples, and are 
therefore strongly related to each other. In the following 
we will consider only the FOT line in the phase diagram. 

M(T) data provide additional information about the vor- 
tex behavior. In the reversible region above Ti rr a clear 
diamagnetic signal is present up to temperatures larger 
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FIG. 2: a) Real and b) imaginary part of the ac-susceptibility 
X(T) of UD-29K in an external field of 3 T. T FO t is deter- 
mined by the peak position in x"(T), which corresponds to 
the maximum slope in x'{T). c) Magnetization data mea- 
sured at Hdc=3 T, after subtraction of a linear background 
taken in the normal state. Below T irr the FC and ZFC M(T) 
curves separate, whereas at Tfot a jump is observed. T C 2 
is estimated by extrapolation (see text). Tfi uct is defined as 
the temperature where the data deviates from the horizontal 
normal state line. 



than T c . This region is characterized by strong fluctu- 
ations and there is no well defined upper critical tem- 
perature T c2 . The temperature Tfi uct , at which diamag- 
netic (superconducting) fluctuations appear, has been de- 
fined as the temperature where the data begin to deviate 
from the horizontal normal state line (see Fig|2t). The 
simplest way to estimate T C 2 is to use the extrapolation 
method based on the linear Abrikosov formula^. The 
transition temperature T c i is derived from the intersec- 
tion of a linear fit with the normal-state horizontal line, 
as shown in Fig[2t. It was shown that this procedure is 
not totally correct for HTSC, where the Abrikosov linear 
region is limited to a small temperature range because of 
the rounding close to T c ^&2&. Indeed in the underdoped 
regime, where fluctuations are larger, using extrapola- 
tion we get unphysical values for the upper critical field 
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FIG. 3: ZFC isothermic magnetization curves for a) OD-31K 
and b) UD-29K. For OD-31K, H p and H sp have been deter- 
mined as indicated by the arrows. For UD-29K, only H sp 
could be observed. The insets show full hysteresis loops with 

Hirr • 



(positive slope of H C 2(T), see Sec. II V|) . However, treat- 
ing the data as proposed by Landau and Ot^S one gets 
more reasonable upper critical lines for all the samples 
(see FigHJ). 

We also performed isothermic ZFC M(H) measurements 
at different temperatures (see Fig|3J|. In the OD sam- 
ples we could observe two peaks in the M(H) curves 
(see FigEK for OD-31K). The first minima H p in the 
OD samples is known to be related to surface^ and/or 
geometrical^ barriers. Due to these barriers the field 
doesn't penetrate the bulk at the lower critical field H c \ 
but only at an higher field H p . The second (and largest) 
minima H sp (second peak) is related to some flux-pinning 
mechanism, although its origin is controversiaiL&iLiil. In 
UD samples only one peak could be observed. We argue 
that this is actually the second peak H sp . The penetra- 
tion field H p is most probably hidden, due to the low 
value of H sp . This interpretation is supported by the 
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in the insets of FigEl The ascending and the descending 
branches meet at Hi rr , whose values are consistent with 
those obtained by FC-ZFC M(T) curves. 
In order to facilitate the analysis and discussion of the 
experimental results, the characteristic fields (H C 2(T), 
Hfi U ct(T), Hfot{T) and H sp (T)) of the four samples 
have been plotted in the H vs T phase diagrams shown 
in Fig^J The magnetic phase diagram of LSCO is usually 
divided in four main phases: 

1. Above the upper critical field H C 2(T), LSCO is in 
the non-superconducting state and the magnetic 
flux is free to enter the crystal homogeneously. 

2. Between H c2 (T) and H FOT {T) (H trr (T)) the mag- 
netic flux is partially expelled from the supercon- 
ductor. The magnetic field is present in the sample 
in the form of vortices which are in a reversible 
regime. In this region the vortices are thermally 
activated and highly dynamic. 

3. Below Hfot{T) (Uj rr (T)) the vortices are in an 
irreversible regime, as can be seen by the difference 
in the FC/ZFC data or in the hysteresis loops. Here 
the vortices are frozen in a lattice (VL), which can 
be directly observed in SANS experimentsiSii&ii. 

4. Below H P (T) (ideally H cl (T)) the system is in the 
Meissner state and the flux is completely expelled 
from the bulk of the sample. 

Indeed we can roughly understand our results in LSCO 
within this description, even though we have some addi- 
tional lines in the phase diagram (e.g. H sp and Hfi uct ). 
The first observation is that the magnetic phase diagrams 
of OD and UD LSCO are qualitatively similar but quanti- 
tatively very different. In particular for the UD samples 
the reversible region is much larger than for OD ones, 
whereas the second peak line occurs at much lower fields. 



FIG. 4: Magnetic phase diagram of the four LSCO samples 
(OD-31K, OD-36K, UD-29K, UD-19K) showing the temper- 
ature dependencies of the second peak field H sp (T), the FOT 
line Hfot(T), the upper critical field H C 2(T) (determined 
by extrapolation and by the scaling procedure), and the field 
Hfi uc t(T) where diamagnetic fluctuations set in. In a)-d) the 
second peak line has been fitted by the power law (Eq.JSJ), 
whereas the FOT line has been fitted by the sublimation 
model (Eq.l|3Jl). In b) we have attempted to fit the second 
peak line by the decoupling theory (Eq.@), while in c) a fit 
of the FOT to the melting theory (Eq.©) is also shown. 



fact that even in the OD samples it is difficult to identify 
Hp at high temperatures close to T c (where H sp occurs 
at low fields). Moreover, very accurate SQUID measure- 
ments on UD-29K clearly showed the presence of two 
minima at H p and H sp even in the underdoped regime 1 -. 
We also performed some full hysteresis loops, as shown 



IV. DISCUSSION 

Before discussing the possible reasons for this strong 
doping dependence of the phase diagram we want to have 
a detailed look at the single lines. 

We start from the upper critical line H C 2(T), which is 
not well defined since fluctuations are very strong near 
T C 2- This is more pronounced in the underdoped regime, 
where diamagnetic fluctuations are present even at tem- 
peratures Tfi uc t much larger than T c . This anomalous be- 
havior in the underdoped regime has also been observed 
in Nernsli 3 ^ 3 ^^ and scanning SQUID microscop}^ ex- 
periments and has been interpreted as being due to 
vortex-like excitations in the pseudogap region. As a con- 
sequence, H C 2(T) as determined by extrapolation has an 
unphysical positive slope. In order to get more reasonable 
upper critical field lines, we used the Landau-Ott scaling 
procedure for magnetization data 30 , taking the values of 
H c2 (0 K) listed in Table [Q The resulting H c2 (T) lines 
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are plotted in Fig0] 

We turn now to the FOT line H F qt(T), which is usually 
identified with the melting lin o 3 i 4 ^ , that is the transition 
of the vortex-solid into a vortex-liquid, in which the VL 
loses its shear modulus. The temperature dependence of 
Hfot(T) is predicted by the melting theory to b o 3 ' 4 '^ : 



H me it{T) — H„ 



T 

1 



(1) 



The prefactor is known to depend almost only on the 
anisotropy of the system. In fact, considering H m ~ 
7~ 2 Tjr 2 A a ^ 4 (A a (, is the in-plane penetration depth) and 
the fact that T~ 2 X~^ is almost constants^, one obtains 
H m ~ 7~ 2 - Fitting our data by this model is not sat- 
isfactory, since we obtain a huge doping dependence of 
the exponent m and the prefactor H m doesn't follow the 
expected 7 dependence (see Table QJ. Moreover, in all 
SANS experiments on HTSC 13 i 40 i 41 the ring-like inten- 
sity expected between Hpor and H c2 for a liquid of 
straight vortices 4 ^ has never been observed. A more 
precise melting theory, still based on the Lindemann 
criterion 4 ^, predicts a more complicated temperature de- 
pendence of the melting line^: 



H me it{T) = 



4c 4 # c2 (0)§(2i 



l + A /l + 4c|f(f -1) 



(2) 



where G = K ( 4ff / Mo ]gg(o)&(o) ) 2 is the Ginzbur S num " 
ber (/i is the permeability of free space, kg is the Boltz- 
mann's constant, H c is the thermodynamic critical field, 
and £ a ;, is the in-plane coherence length), B«5.6 and Cj, 
is the Lindemann number. However, even this formula 
doesn't describe our data very well, since the fitted curves 
are unsatisfactory (see for example FigQJi for UD-29K), 
c L is doping dependent and in some cases higher than 
the expected values (c L ~0.1-0.2). 

An alternative model to the melting transition is given by 
the sublimation theory^, based on the strong anisotropy 
originating from the layered structure intrinsic to all 
HTSC. Within this scenario the melting is accompanied 
by the simultaneous decoupling of the vortex lines into 
2D pancake vortices (vortex gas). The phenomenological 
scaling law which applies to all HTSC is given by: 



H sm {T)[Oe] 



2.85 ■7 _2 ,s- 1 



Tc 
T 



1 



(3) 



where s is the distance between the Cu02 layers 
(6.6T0 -8 cm in LSCO). This formula has been used in 
order to explain the FOT transition and nicely fits our 
data. 7 is the only free parameter, and the fitted values 
are in good agreement with the measured values of the 
anisotropy (see Table QJ. 

It remains to discuss the second peak line H sp (T) which 
has been explained on the basis of the thermal decou- 
pling theoryi^Siii, which predicts the suppression of 



TABLE I: Characteristic parameters for LSCO as a function 
of the Sr concentration x. The values of the upper critical 
field H C 2(0 K)-, of the penetration depth A a »", and of the 



anisotropy 



have been extrapolated from experimental 



values found in the literature. H m and m have been obtained 
by fitting the data by Eq.Q, the Lindemann number cl using 
Eq.@. 7su6i and 7d ec are the anisotropics obtained by fitting 
our data using the sublimation, respectively decoupling mod- 
els. Finally, the exponent n has been obtained by fitting the 
data with Eq.©. 
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long-range order in the direction of the applied field due 
to thermal fluctuations. The expected temperature de- 
pendence is 4 ^ 



ruA T i n ( ^ - 1 



(4) 



with H* = $g/(167r 3 efc B /u S7 2 T c A a6 (0) 2 ), where <J> is 
the flux quantum and e ~ 2.718 is the exponential num- 
ber. This function doesn't fit well our data, as shown in 
FigEI 3 f° r OD-36K. Moreover, the estimated values for 
7, obtained by substituting the known values of s, T c 
and A a fc(0) in the theoretical expression for B* , are not 
satisfactory compared to the experimental values (see Ta- 
ble P) . Moreover recent SANS measurements^ 3 - indicate 
that the diffraction signal from the vortex lattice persists 
up to Hfot(T) and therefore discredit the decoupling 
theory. The origin of the second peak is most proba- 
bly to be found in some change of the pinning mech- 
anism. It has been often suggested that this feature is 
related to the transition to a more disordered vortex glass 
phas o 9 ' 10 , and very recent experimental results confirm 
this interpretationAi* 1 ^. Our experimental data are bet- 
ter fitted by a power law 4 ^ 



H sp (T) = H ■ 1 



T 

TV. 



(5) 



as can be seen in Fig 0| and Fig[5j The value of the expo- 
nent is close to n=2 in all samples (see FigEland Table 
P). Interestingly the value of H seems to be propor- 
tional to 7~ 3 , even though (up to our knowledge) no 
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FIG. 5: Temperature dependence of H sp ■ 7 3 plotted in a 
double logarithmic scale. All the data measured in samples 
with different doping levels collapse on one line with slope ~ 2. 
This indicates that the power law (Eq.l0) has an exponent 
n ~ 2 and Ho oc 7~ 3 . We have used the values of 7 obtained 
by fitting our data using Eq.@ (jsuu in Table QJ. 



theory predicts such a 7 dependence. However, a large 
anisotropy naturally renders the vortex system more sus- 
ceptible to disorder. The observed anisotropy depen- 



dence of H sp (T) is therefore in qualitative agreement 
with a scenario where the second peak line is related to 
a field-induced vortex glass transition. 



V. CONCLUSION 

A first look at the magnetic phase diagrams shown in 
FigEl could indicate that the vortex matter in LSCO is 
strongly doping dependent. This is true from a quanti- 
tative point of view, but qualitatively all samples display 
the same transitions (second peak, irreversibility, FOT 
and upper critical lines). The quantitative doping de- 
pendence of the magnetic phase diagram can mainly be 
explained by the different degree of anisotropy: H sp is 
found to be proportional to 7~ 3 and Hfot to 7~ 2 . The 
interpretation of the second peak in LSCO is still con- 
troversial but our data seem to favor the vortex glass 
scenario, whereas the FOT line is consistent to the subli- 
mation theory rather than to the melting theory. More- 
over, strong superconducting fluctuations above T c have 
been observe in the underdoped regime. 
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